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Preface 

This  report  is  a  stiinmary  of  activities  carried  on  at  the 
Institute  of  Mathematical  Sciences  at  New  York  University  from 
July  1957  to  July  19580   It  contains  some  personnel  data,  a 
summary  of  scientific  progress  during  the  past  year,  a  list  of 
reports  and  papers  Issued,  and  finally  a  list  of  references. 
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Personnel 

The  average  number  of  equivalent  full-time  scientific  staff 
members  is  16.6  theoreticians  and  k   experimenters.  The  number  of 
individuals  involved  is  somewhat  greater.   Of  these,  8  have  facul- 
ty rank  or  equivalent  senior  positions.  There  are  12  Ph.D. 's  in 
mathematics,  and  6  with  Ph.D. 's  in  physics  or  other  fields.  Also 
employed  are  8  research  assistants.  During  the  year  5  consultants 
have  been  used. 


Shock  Theory  [l8,26be,27ab] 

The  evidence  has  become  increasingly  strong  that  there  do 
exist  steady  shock  profiles  with  an  irreversible  transition  from 
one  constant  state  to  another  without  the  intervention  of  collisions 

as  a  damping  mechanism.  Most  of  the  investigation  has  concerned  a 

* 
range  of  parameters  in  which  a  low  temperature,  low  p  plasma  is 

heated  to  a  much  higher  temperature  at  a  moderate  p  by  a  shock 
of  moderate  strength.    It  was  established  previously  that  such  a 
shock,  provided  that  it  is  observable  within  reasonable  apparatus 
dimensions,  heats  the  ions  almost  exclusively.  The  wave  length  of 
the  oscillating  shock  profile  is  on  the  order  of  'K=^m/\me           and 
the  entire  thickness  of  the  transition  can  range  from  several  times  "K 
to  many  times  X,  depending  on  the  parameters.   A  considerable  amount 
of  analytical  work  has  been  done  involving  all  the  general  features. 
Exact  quantitative  results  require  numerical  computation  of  an  ex- 
tremely intricate  and  non-standard  type.  This  program  is  xinderway. 
— _ _ — — __ — _ — _ _ _ _ 

A  representative  value  of  X  is  .02  cm.  at  n=10  vcm  » 


«»* 
«« 


Strength  is  measured  by  magnetic  field  ratio, 
P=2m.p/B  =  gas  energy /magnetic  energy, 
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Runaway  Ions  and  Thermonuclear  Reaction  Rates  [12,19,26d] 

It  was  remarked  at  the  Princeton  meeting  (Oct. '57)  that, 
because  of  a  dip  In  the  Ion  dynamical  friction  curve  occurring  at 
fast  Ion  speeds  Intermediate  between  the  mean  Ion  and  mean  elec- 
tron speed,  one  might  expect  semi -runaway  Ions  to  develop  on 
application  of  a  moderate  electric  field.   Subsequent  quantitative 
analysis  showed*-   '^*    ■'  that  no  such  effect  exists,  because 
the  Jon  distribution  remains  very  accurately  Maxwellian  even  with 
an  applied  E,  and  the  dip  in  the  friction  curve  remains  far  out 
in  the  tail  of  the  ion  distribution »   However,  computation  of  the 
(small)  perturbation  of  the  ion  distribution  from  Maxwellian  leads 
to  a  corrected  thermonuclear  reaction  rate.  The  Increased  rate 
can  be  quite  large  under  certain  circumstances,  and  a  center  of 
mass  motion  opposite  to  the  applied  E  is  predicted.   It  is  not  yet 
clear  whether  these  effects  are  observable.  From  this  analysis, 
it  appears  that  actual  runaway  ions  (with  speeds  greater  than  the 
mean  electron  speed)  cannot  appear  under  any  conditions  which  are 
experimentally  realizable.   In  computing  the  perturbation  of  the 
ion  distribution,  mathematical  methods  of  computing  the  tail  end 
of  the  distrubutlon  were  developed  which  may  have  many  other 
applications  (e.g.  to  co  mpu';ation  of  runaway  electrons). 

Theory  of  the  Bojtzmann  Equation 

A  very  complete  set  of  transport  coefficients  has  been  com- 
puted for  a  fully  ionized  gas.^^'  ^*-^  •'  These  include  many  more 
coefficients  than  those  of  the  Chapman-Ens -cog  theory,  e.g.  a  2  x  2 
matrix  of  heat  conductivities  (connecting  electron  and  ion  tempera- 
tures with  electron  and  ion  carried  heat  flows),  a  similar  multi- 

-  k  - 


pllclty  of  viscosity  coefficients,  thermal  diffusion  coefficients, 
together  with  all  the  additional  "Hall"  effects  and  the  distinc- 
tion between  parallel  and  perpendicular  values  of  each  coefficient 
in  the  presence  of  a  magnetic  field.  These  coefficients  include 
all  the  standard  values  of  electrical  conductivity.  Much  work  has 
to  be  done  to  assess  these  results  and  discover  under  what  circiim- 
stances  the  Chapman-Enskog  -"^alues  are  inapplicable. 

From  the  moment  analysis  of  the  Boltzmann  equation  one  can 
derive  a  set  of  macroscopic  fluid-like  equations  which  contains 
the  effect  of  finite  Larm  or  radius  together  with  anis  tropic 
stress  tensor  and  heat  flow  and  current  flow  along  the  field  lines. 
These  equations  should  be  considerably  more  accurate  than  the 
guiding  center  (or  Chev;  — Goldberger-Low)  equations.  They  are 
probably  easier  to  apply  than  the  zero  Larmor  radius  self-consis- 
tent field  technique  and,  in  addition,  take  into  accotint  some  of 
the  effects  of  finite  Larmor  radius.   So  far,  they  have  been 
applied  to  only  the  shock  problem  with,  as  yet,  inconclusive 
results^ ^^^. 

The  validity  of  the  Boltzmann  and  Fokker-Planck  equations 
has  been  studied '• -^"^ ^ .   It  appears  that,  despite  the  fact  that  the 
two-particle  correlation  distance  (viz.  the  Debye  length)  is  much 
larger  than  the  mean  interparticle  separation  distance,  the  one- 
particle  Boltzmann  equation  can  be  used  (raising  one's-eelf  by 
one ' s  boot  straps )  to  compute  the  correct  non-equilibrium  shielding 
potential.   Preliminary  results  have  been  computed.   A  tentative 

IF ~— — -~ ^ ~~~ ^ ~~ ~    ~~      ' 

Preliminary  results  were  given  by  M.  H.  Rose  in  Gatlinberg, 
J\ane  I956. 
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correction  of  the  conventional  Boltzmann  equation  (with  Coulomb 
potential  cut  off  at  the  Debye  distance)  using  the  spherically 
symmetric  shielding  potential,  e'yr,  has  been  computed'-  ■^' 

Spiraling  of  Particles  [11] 

The  constancy  of  the  magnetic  moment  p.  was  demonstrated 
formally  by  Alfven  to  first  order,  by  Hellwig  to  second  order,  and 
by  Kruskal  to  all  orders.  These  formal  expansions  have  been 
shown  to  be  actually  convergent  to  solutions  of  the  equations  of 
motion.  The  problem  is  to  estimate  the  error  made  in  terminating 
the  expansion  at  a  given  point  and  is  a  difficult  one  because  of 
the  possibility  of  cummulation  of  error  after  many  Larmor  cycles. 

Self  Consistent  Field  Problems 

Techniques  have  been  developed  for  solving  by  approximate 
means  problems  of  much  more  generality  than  are  accessible  by 
exact  methods.   A  specific  problem  that  has  been  treated  this  way 
is  the  rotating  mirror  or  homopolar  machine.   Once  the  self- 
consistent  equilibrium  solution  has  been  found,  the  particle 
losses  can  be  computed.   Some  results  were  presented  at  the 
theoretical  meeting  at  N.Y.U.in  May'-^'''®^  The  techniques  involve 
a  combination  of  the  pre\lously  developed  generalized  adiabatic 
invariance  theory^ -^-^ ^ ,  the  guiding  center  fluid  theory "-^  \   to- 
gether with  certain  other  techniques  for  treating  rotating 

r2Ql 
systems'-  ^\ 

Stability 

A  large  ntamber  of  mirror  machine  geometries  with  end  plates 

[17 
were  analyzed  for  stability  using  a  simple  su  f ficient  condition'  ' 

26c, 27d]^   Stable  configurations  were  found  with  interesting  values 
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of  p.  The  most  interesting  ranges  of  the  parameters  (viz.  mirror 
ratio  and  aspect  ratio)  will  be  recomputed  using  the  more  compli- 
cated necessary  and  sufficient  analysis. 

The  fact  that  a  mirror  machine  has  open  ends,  past  which 
waves  can  pass  remaining  energy  from  the  system,  is  a  stabilizing 
influence  that  has  not  been  previously  analyzed.   A  rough  analysis 
of  this  problem  yields  9  finite  amplitude  of  oscillation  in  sys- 
tems which  would  otherwise  be  called  unstable. 

Cusped  Geometries 

Methods  of  computing  three  dimensional  cusped  geometries  both 
analytically  and  numerically  are  being  developed.   Particle  losses 
from  low  p  picket  fence  geometries  are  being  computed  in  order 
to  have  precise  results  for  comparison  with  cusped  geometries. 
In  many  respects,  the  theory  has  been  carried  as  far  as  that  which 
is  available  for  other  devices,  and  no  large  theoretical  program 
is  warranted  without  experimental  guidance. 

Experimental  Program 

The  laboratory  was  transferred  to  new  quarters  with  greater 
space. 

An  automatic  charging  and  firing  system  was  designed  and  built 
with  improved  safety  features,  permitting  experiments  to  be  per- 
formed in  five  seconds  instead  of  the  previous  100  seconds  or  more. 
This  has  materially  increased  capacitor  lifetime.   In  addition  , 
the  main  capacitor  bank  was  redesigned  for  horizontal  operation  and 
has  been  fired  repeatedly  at  three  kilovolts  (maximum  voltage)  with- 
out breakdown. 
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A  new  vacuvim  system  was  put  into  operation,  permitting 
experiment  pressures  down  to  10~  mm  Hg.   A  new  mass  spectro- 
meter leak  detector  system  has  been  installed  also. 

New  rail  and  coaxial  sources  have  been  constructed,  employed 
up  to  ten  kilovolts;  the  coaxial  source  has  evidenced  more  uni- 
form plasma  burst  alignment. 

With  this  equipment,  successful  plasma  constrictions  (pipes) 
up  to  one  meter  have  been  obtained  in  both  two  inch  and  one  inch 
straight  glass  tubes.   Conversely,  U-bend  studies  show  that  in 
regions  of  high  field  c\irvature  the  plasma  flow  curvature  remains 
small.   However,  a  tube  of  larger  radius  of  curvature  has  been 
constructed  which  will  employ  greatly  increased  fields. 

Time-of=flight  measurements  were  performed  on  the  plasma  in 
the  one  inch  tube.   Results  indicated  plasma  velocities  of  over 
two  centimeters  per  microsecond,  essentially  supersonic. 

A  large  series  of  experiments  were  perfonned  in  the  one 
inch  tube,  shooting  plasma  bursts  along  axial  fields,  up  to  twenty 
kilogauss,  which  has  been  distorted  or  diverted  by  brass  obstacles. 
The  plasma  appeared  to  follow  the  diverted  field  lines  as  predicted 
by  theory. 

A  report''   ■'  was  prepared  and  distributed  of  all  of  the  im- 
portant results  obtained  since  July  1,  19^1. 

A  new  sixty  kilovolt  coaxial  plasma  source  (cannon)  was 
recently  completed,  together  with  its  own  vacuum  spark  gap  trigger 
assembly.   For  the  present,  it  is  being  employed  at  twelve  kilovolts 
(twenty  times  the  previous  source  energy)  to  gather  preliminary 
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data,  and  to  avoid  the  possibility  of  destructive  testing. 

A  new  two  inch  straight  experiment  tube  was  constructed  on 
which  the  center  third  of  the  main  field  solenoid  was  wound  at 
four  times  the  pitch  of  the  outer  thirds,  thus  (theoretically) 
increasing  the  fl\ix  density  in  this  region  four-fold.  This  tube 
was  assembled  with  the  new  cannon  source.   Preliminary  results 
show  conclusively  that  the  plasma  converges  in  the  central  region 
where  the  field  lines  would  be  expected  to  converge. 

A  detailed  tirae-of-flight  study  has  been  almost  completed  on 
the  above  system.   Preliminary  results  suggest  that  the  plasma  burst 
expands  longitudinally  since  the  leading  and  trailing  edges  evidence 
different  velocities.   If  this  velocity  difference  were  interpreted 
as  temperature  it  would  Indicate  thermal  energies  in  excess  of  ten, 
and  possibly  one  hundred,  volts. 
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Reports  and  Papers  Issued,  July  1957  -  July  1958 

MH-V  (Ref.  [5]) 

This  report  contains  material  which  was  to  a  large  extent 
carried  out  in  195^  and  summarized  at  the  Princeton  meeting, 
1954.  The  mathematical  theory  of  Maxwell's 

equations  without  displacement  current  is  developed  extensively 
and  is  found  to  be  considerably  simpler  and  more  easily  applied 
than  that  of  the  full  Maxwell's  equations,  e.g.  in  the  transient 
theory  of  the  pinch.   This  theory  (which  is  Galilean  invariant)  is, 
in  a  sense,  more  appropriate  thati    the  theory  of  the  full 
equations,  when  used  together  with  non-relativistic  fluid  equations. 

MH-VI  (Ref.  [6]) 

This  report  contains  background  material  in  magnetohydro- 
dynamics  derived  and  collected  for  easy  reference  in  later  reports. 

MH-VII  (Ref.  [7]) 

In  this  paper  a  number  of  exact  mathematical  analogues  of 
magnetohydrodynamic  problems  with  conventional  fluid  dynamical 
problems  are  presented.  This  makes  available  a  wealth  of  litera- 
ture and  experience  for  use  in  these  special  problems.   This 
material  was  developed,  in  large  part,  in  195^ >  and  was  Issued  in 
a  mimeographed  set  of  NYU  lecture  notes  in  195^' 

MH-VIII  (Ref.  [8]) 

This  contains  an  elaboration  and  extension  of  the  material 
contained  in  LAMS-2105,  issued  in  1957.  The  theory  of  the  com- 
plete set  of  nonlinear  classical  magnetohydrodynamic  equations 
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(including  Maxwell's  equations  without  displacement  current)  is 
developed  and  many  techniques  of  gas  dynamics  are  found  to  be 
directly  applicable.   Small  disturbances  and  hydroraagnetic  shocks 
are  treated,  as  well  as  the  theory  of  "simple  waves". 

Ref.  [11] 

Discussed  above. 

Ref.  [12] 

Discussed  above.  There  is  also  included  an  analysis  of  the 
relation  between  the  Boltzmann  and  FoKker-Planck  equations. 

Ref.  [13] 

Discussed  above. 

Ref.  [14] 

The  magnetohydrodynamic  theory  of  the  full  Maxwell's  equations 
coupled  with  a  set  of  relativistic  fluid  equations  is  shown  to  be 
mathematically  equivalent  to  the  theory  of  the  conventional 
equations  (MH-VIII). 

Ref.  [15] 

Discussed  above. 

Ref.  [16] 

Discussed  above. 


-  11 


Papers  for  the  Second  U.  N.  International  Conference  on  the 
Peaceful  Uses  of  Atomic  Energy  In  Geneva,  Sept.  1958 


Ref.  [17] 

Contains  a  complete  acco\int  of  the  theory  of  stability  of  free 
boundary  equilibria  treated  by  variational  methods  [35*27d,246e] , 
also  the  analysis  of  a  special  problem  with  rotating  magnetic 
fields;  also  the  mirror  machine  stability  computation  referred  to 
above. 

Ref.  [18] 

Discussed  above. 

Ref.  [19] 

Discussed  above. 

Ref.  [20] 

Contains  a  summary  of  the  theory  of  magnetohydrodynamic  equi- 
librium (of  which  force  free  fields  are  a  special  case).  The 
results  presented  embrace  considerably  more  material  than  exists 
in  the  previously  published  literature. 

Ref.  [21] 

Contains  a  condensation  of  the  theory  of  cusped  geometries. 
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